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Q ; Abstract 

A systematic study of the temperature dependence of the shapes and pair- 
ing gaps of some isotopes in the rare-earth region is made in the relativistic 

o 

■ Hartree-BCS theory. Thermal response to these nuclei is always found to lead 

|Q ■ to a phase transition from the superfluid to the normal phase at a temper- 

ature Ta ~ 0.4 — 0.8 MeV and a shape transition from prolate to spherical 
q ■ shapes at T c ~ 1.0 — 2.5 MeV. These shape transition temperatures are appre- 

ciably higher than the corresponding ones calculated in the non-relativistic 
: ^ ■ framework with the pairing plus quadrupole interaction. Study of nuclei with 

X 

continued addition of neutron pairs for a given isotope shows that with in- 
creased ground state deformation, the transition to the spherical shape is de- 
layed in temperature. A strong linear correlation between Ta and the ground 
state pairing gap A is observed; a well-marked linear correlation between T c 
and the ground state quadrupole defromation 02 is also seen. The thermal 
evolution of the hexadecapole deformation is further presented in the paper. 
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I. INTRODUCTION 



Heating can have a profound effect on nuclear shapes, causing a variety of shape transi- 
tions. Experimentally, such responses to the thermal excitations have been studied from the 
shapes of the giant dipole resonances(GDR) built on excited states Theoretically, they 

have been studied earlier in a finite temperature non-relativistic microscopic Hartree-Fock 
and Hartree-Fock Bogoliubov (HFB) framework with a pairing plus quadrupole 
(P+Q) interaction. For the nuclei studied, it has been found that while the superfluid nu- 
clear phase has transition to the normal phase typically at T~ 0.5 MeV, the deformed shapes 
have transition to the spherical ones at a higher temperature, mostly between T~ 1.0 to 1.8 
MeV for rare-earth nuclei. These calculations, however, have some limitations; they employ 
a simplistic model Hamiltonian in a limited model space, an inert core is assumed, moreover, 
the Coulomb interaction has not been taken into account realistically. The understanding 
of the universal pattern of the mean field shape evolution with temperature has also been 
tried in a macroscopic approach ||iO| , |TTf commonly referred to as the Landau theory of phase 
transition. A quantitative estimate of the persistence of the ground state deformation |12| 
with temperature is however seen to be missing in some cases. 



Recently, we have undertaken a study |T3[ of the thermal evolution of nuclear properties, 
particularly the phase transition in the nuclear shape and the superfluidity in the relativistic 
mean field (RMF) theory. The pairing effects have been included in the BCS approximation. 
The RMF theory [|T^-[T6[| has proved to be an extremely powerful tool in explaining the gross 
properties of nuclei over the entire periodic table. In contrast to the non-relativistic models, 
this theory employs a single set of parameters to explain all these properties. Moreover, 
in such calculations , the model space used is sufficiently large and all the nucleons are 
treated on equal footing. The calculations reported in Ref. [H| are performed for only two 
rare-earth nuclei, namely, im Er and 170 Er. It is found that the phase transition for the 
nuclear shape from the prolate to the spherical occurs at a temperature significantly higher 
than that obtained in the (P+Q) model and that the transition is relatively smooth. A very 
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recent calculation by Egido et al |17J in a nonrelativistic approach but with the realistic 
Gogny force reveals that the characteristics of the nuclear shape transition are very similar 
to those obtained in the RMF theory. We have therefore undertaken a more systematic 
study of the shape transition for the rare-earth nuclei in the relativistic mean field theory in 
the present paper. For this purpose, we have considered various even-even isotopes of Sm, 
Gd and Dy. A more quantitative study, particularly of the shape transitions, calls for the 
inclusion of thermal fluctuations |13| , [T8|JT9| . However, this is too computer intensive and are 
not included in the present paper. 

The organization of the paper is as follows: we discuss the theoretical framework briefly 
in section II. The results and discussions are presented in section III and the concluding 
remarks are given in section IV. 

II. FORMALISM 



We employ the nonlinear o — uj — p version of the RMF theory ||15|| . The Lagrangian 
density for the nucleon-meson many body system is taken as 

£ = (jyfy -M)% + -d^ad^a - U(a) - g a %a% 

-gjtltfp^i - -F^F^ - e^il^pU^. (1) 
The arrows indicate isovector quantities. The mesons included in the description are the 



a, uj and p mesons. For an appropriate description of the nuclear surfaces ||20|| , a non-linear 



scalar self-interaction term U(a) of the a meson is included in the Lagrangian 

U(a) = \mla* + ^ + \g^ A . (2) 

The meson masses are given by m a , m w and m p , the nucleon mass is M and g a , g w , g p 
and e 2 /47r = 1/137 are the coupling constants for the mesons and the photon. The field 
tensors for the vector mesons uj and p are given by Vt^ v and R^, for the electromagnetic 



field, it is F^ u . Recourse to variational principle followed by the mean field approximation 
treating the fields as c— numbers results in the Dirac equation for the nucleon and the Klein- 
Gordon type equations for the mesons and the photon. For the static case, along with the 
time-reversal invariance and charge conservation, the equations get simplified. The resulting 
equations, known as RMF equations, have the following form. The Dirac equation for the 
nucleon is 

{-ia • V + V(r) +{3[M + 5(r)]}^ = e^, (3) 
where V(r) represents the vector potential 

V(r) = ^o(r) + g p T 3 p (v) + e iI_^U (r), (4) 
and S(r) is the scalar potential 

S(r)=g a a(r), (5) 

which contributes to the effective mass as 

M*(r) =M + S(t). (6) 

The Klein-Gordon equations for the mesons and the electromagnetic fields with the 
nucleon densities as sources are 

{-A + ml}a( r ) = -g aPs (r) - g 2 a 2 (r) - g 3 a 3 (r) , (7) 

{-A + m ^ (r) = ^(r), (8) 

{-A + m 2 p }p (r)= g p p 3 (r), (9) 

-AA)(r) = ep c (r). (10) 
The corresponding densities are 

p s = Y^m^^i, 

i 
i 
i 

p c = £^!^^*, (ii) 



Here the sums are taken over the particle states only, i.e., the negative-energy states are 
neglected. The partial occupancy (ni) at finite temperature in the BCS approximation is 



1 

Hi = - 

2 



l-^(l-2/(^T)) 



(12) 



with f(ii,T) = 1/(1 + e €4//T ); = \J (e« — A) 2 + A 2 is the quasiparticle energy where q is 
the single-particle energy for the state i. The chemical potential A for protons (neutrons) is 
obtained from the requirement 



5> = Z(JV) 



(13) 



The sum is taken over proton (neutron) states. The gap parameter A is obtained by min- 
imising the free energy 



F = E-TS, 



(14) 



where 



E(T) = J2 Urk + E a + E aNL + E UJ + E p + E c + E pair + E c , m , - AM, 



(15) 



and 



S = - y Z\filnfi + (l-fi)ln(l-f i )], 



(16) 



with 
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\g* j d 3 rp s {r)(r{r), 
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(17) 
(18) 
(19) 
(20) 
(21) 
(22) 
(23) 



Here G and A are the pairing strength and the mass number respectively. The single-particle 
energies and the fields appearing in eqs. (|I5D - (^1|) are obtained from the self-consistent 
solution of eqs. (||) - (|T0|). The temperature dependent occupancies of the fermions induce 
temperature dependence in the bosonic fields through the source terms as seen from eqs. 

m-m- 



We generate these self-consistent solutions using the basis expansion method |15| , |21| ; 
this yields the quadrupole deformation j3 2 , the hexadecapole deformation /3 4 and the proton 
and neutron pairing gaps A p and A n as a function of temperature. 

III. RESULTS AND DISCUSSIONS 

We have chosen even-even isotopes of the nuclei Sm, Gd and Dy for the study of the 
pairing and shape transitions. The NL3 parameter set is chosen for the values of the coupling 
constants and the masses of the mesons and the nucleons. This parameter set reproduces best 
the ground state as well as the compression properties of finite nuclei simultaneously |22f ; 
however, it has already been reported in Ref. |13|] that the results for shape transition are not 
that sensitive to the choice of the parameter set. The pairing gaps A° and A° for neutrons 
and protons for a nucleus in the ground state are determined from the experimental odd- 



even mass differences [23]. The single-particle states are calculated using spherical oscillator 
basis with twelve shells. The values of the chemical potential and the pairing gap at a given 
temperature are determined using all the single particle states upto 2huo (the model space) 
above the Fermi surface without assuming any core. 

At finite temperature, because of the partial occupancies of nucleons above the Fermi 
surface, it is in principle necessary to have a larger basis and an extended model space. It 



is further necessary to take effects due to continuum into account In order to check 

the convergence of the calculations, we have enlarged the basis space from twelve shells to 
twenty shells and have extended the model space to include single-particle states upto 3Hujq 
above the Fermi surface. For this extended model space, the pairing strength G is adjusted 
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to reproduce the ground state pairing gap. The changes in the values of the observables 
are found to be insignificant due to this extension of the basis and model space even at 
the highest temperature of our interest (T ~ 3.0 MeV). To estimate the importance of the 
continuum corrections on the observables we report here, we calculated the occupancy 
of the single particle states with positive energy. For T < 1 MeV, practically there is no 
particle in the positive energy states = 0) and at the highest temperature of interest 
studied here (T = 2.7 MeV), n {+) / A = 0.011 which is very small. It is therefore expected 
that continuum corrections may not play an important role in the temperature range we 



study. Calculations of nuclear level density in earlier studies ||T3|j25[1 , have shown that the 
continuum corrections are not important for T upto ~ 3 MeV. The continuum effects may 
grow stronger for T > 3 MeV, however, this is beyond the shape transition temperatures 
and so we have not taken this into account. 

The temperature evolution of the quadrupole deformation /?2 with neutron number N = 
86 and 88 for the systems Sm, Gd and Dy are displayed in Fig. 1. It is well known that 
addition of nucleons beyond the closed shell gives nuclei progressive prolate deformation upto 
around the middle of the next shell closure. This is reflected in the figure for all the isotopes 
and isotones. It is also seen that the critical temperature increases with addition of nucleons 
for these systems. It is not immediately apparent whether there is a close correlation between 
the ground state deformation fl\ and the critical temperature T c for shape transition; we 
come back to this issue later. In the top panel of this figure, particularly for 148 Sm, it is 
seen that the deformation increases a little with temperature before finally falling to zero. 
This is due to the delicate balance between the temperature dependence of the pairing 
force and the nuclear interaction as derived from the RMF theory. The dramatic build- 
up of a deformation in a temperature window for this nucleus as seen earlier |||7[] in the 
non-relativistic framework is absent in our calculations. 

The temperature dependence of the pairing gaps A p and A n for protons and neutrons 
for the two isotopes each of Sm, Gd and Dy are shown in Figs. 2 and 3. The pairing gaps 
decrease monotonically with temperature, vanishing at T ~ 0.6 — 0.7 MeV for neutrons and 



at ~ 0.65 — 0.85 MeV for protons. The sudden collapse of the pairing gap and the nuclear 
deformation at some specified temperatures signifies phase transitions; these correspond to 
transition from the superfluid nuclear phase to the normal phase and a transition from the 
deformed shape to the spherical shape, respectively. These transitions show up as bumps 
(displayed in Fig. 4 for 148 Sm and 150 Sm) in the temperature evolution of the heat capacity 
defined as 

C{T) = (24) 

where E* is the excitation energy of the nucleus in question. At a temperature Ta ~ 0.6 
MeV, the twin peaks are seen for both the nuclei referring to the dissolution of the neutron 
and proton pairing gaps. These are the characteristic signatures of second order phase 
transition from the superfluid to the normal phase. A somewhat more prominent bump is 
seen at a temperature T c ~ 1.15 MeV for 148 Sm (upper panel). This corresponds to the 
nuclear shape transition. Addition of two neutrons (lower panel for 150 Sm.) shifts the shape 
transition temperature to T c ~ 1.6 MeV. This is possibly due to the larger ground state 
deformation of the 150 Sm nucleus. 

From the study of the hot 148 Sm and 150 Sm nuclei , it was conjectured earlier |7J] that 
addition of two neutrons might increase the critical temperature for deformation collapse. To 
test this conjecture in detail, we have calculated the ground state quadrupole deformation 
and shape transition temperatures for a host of even-even Sm isotopes. The results are 
displayed in Fig. 5. In the vicinity of the closed shell (N = 82), the ground state quadrupole 
deformation (3% increases fast with the addition of two neutrons as seen from the top panel 
of Fig. 5. As the neutron number approaches the mid shell, the deformation levels off, and 
then, as is well known, switches over to the oblate shape [^6] (not shown in the figure). The 
critical temperature (T c ) for the shape transition also increases with neutron pair addition 
as seen in the bottom panel of Fig. 5. The functional behavior of (3® and T c with neutron 
number N are found to be very similar . Indeed there is a strong correlation between (3% and 
T c as displayed in Fig. 6. The filled circles refer to the results from the present calculation; 
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they can be fairly well fitted with a straight line 

T c = 7.75/3°. (25) 

The fit is obtained using results from the twelve different isotopes of Sm, Gd, Dy and Er 
nuclei. The results calculated J7|||]27J in the (P+Q) model are also presented in the figure 
(open squares) for a comparison with those obtained from the RMF theory. The T c — (3® 
correlation is then also found to be approximately linear with a smaller slope. The shape 
transition temperatures obtained in the present calculations are somewhat higher compared 
to those obtained in the (P + Q) model; one may be inclined to attribute this difference to 
the higher effective mass in the (P + Q) model. However, calculations with different sets of 
field parameters in the RMF theory with considerably different values of the effective masses 
yield conflicting results. We have done calculations for 150 Sm with the parameter sets HS 
and NL2 which yield very different M*/M (0.54 and 0.67), but the T c comes out to be 
1.45 and 1.75 MeV, respectively, contrary to the simple-minded expectations. Therefore no 
simple explanation for the higher values of T c in the RMF theory is obvious. 

How does the pairing transition temperature T A depend on the pairing gap at T = 0? 
Intuitively one would expect the collapse of nuclear superfluidity at a larger temperature 
if the ground state pairing gap A° is larger. To study it quantitatively, we have done 
calculations for the different isotopes of Sm. The results are shown in Fig. 7. In its 
upper panel, we display the experimental neutron and proton ground state pairing gaps as 
a function of neutron number for the Sm isotopes; in the lower panel, the pairing transition 
temperatures T A are displayed. An extremely strong correlation between T A and A° for 
both neutrons and protons is seen; this is manifest in the linear relationship 

T Ap = 0.56A;, T An = 0.60A°, (26) 

for both neutrons and protons which is also shown in Fig. 8. The points in the figure include, 
in addition to Sm, results from Gd, Dy and Er isotopes. The relation between Ta and A° 
is in very close agreement to that obtained in the nonrelativistic (P+Q) model @. 
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The hexadecapole deformation 04, if any, also collapses at the same shape transition 
temperature T c as the quadrupole deformation. In Fig. 9, we plot the hexadecapole moment 
(a measure of f3^) as a function of temperature for 148 SVn and 150 SVn. Addition of neutrons 
gives larger ground state fil- The deformation increases smoothly at low temperature upto 
T ~ 0.7 MeV and then collapses to zero at T c . The initial enhancement of this deformation is 
related to the weakening of the pairing correlations with temperature. Such an enhancement 
is also seen in a finite temperature Hartree-Fock Bogoliubov (FTHFB) calculation using the 



finite range density- dependent Gogny force 17 



IV. CONCLUSIONS 

The relativistic mean field theory has been applied to understand properties of some rare- 
earth even-even nuclei at finite temperature. Pairing effects have been included through the 
BCS approximation. Focus is made mainly on the temperature-induced transition from 
the nuclear superfluid phase to the normal phase and also on the shape transition from a 
deformed shape to a spherical one. To find out the systematics of the dependence of the 
pairing and shape transition temperatures on the values of the ground state pairing gap 
and the ground state deformation, we have done calculations for several isotopes of different 
rare-earth nuclei. We find that there is a linear correlation between the above mentioned 
transition temperatures and the equilibrium values of the pairing gap and deformation at 
zero temperature. The linear relationship is extremely good for the pairing gap and quite 
fair for the deformation. In the range of nuclei that are studied here, it is indeed possible to 
estimated very closely the value of the pairing transition temperature Ta for both neutrons 
and protons, as the ground state pairing gap can be calculated from the systematics of 
binding energy. Similarly, it appears that the shape transition temperatures can be well 
estimated since the ground state quadrupole deformations can be experimentally extracted. 
The transition temperatutres Ta are not too different from those calculated earlier in a non- 
relativistic framework; the shape transition temperatures T c however seem to be higher than 
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the corresponding values calculated in the (P+Q) model. These higher values of the shape 
transition temperatures are however found to be very compatible with the ones obtained 
from the realistic Gogny force. 
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Figure Captions 

Fig. 1 The evolution of the quadrupole deformation j3 2 as a function of temperature for the 
Sm, Gd and Dy isotopes. 

Fig. 2 Temperature evolution of the proton pairing gap for the systems indicated. 

Fig. 3 Same as in Fig. 2, but for the neutron pairing gap. 

Fig. 4 Variation of heat capacity as a function of temperature for 148 Sm and 150 Sm. 

Fig. 5 (a) Variation of the ground state deformation fi\ with neutron number N for the 
isotopes of Sm and (b) the shape transition temperature T c for different isotopes of 
Sm. 

Fig. 6 The shape transition temperature T c plotted as a function of the ground state 
quadrupole deformation fi\ for different rare-earth nuclei. The full circles refer to 
the results from the RMF theory, the full line is the linear fit to these. The open 
squares refer to those from the (P+Q) model and the dashed line is the linear fit to 
these points. 

Fig. 7 (a) The ground state neutron and proton pairing gaps as a function of neutron number 
for the Sm isotopes and (b) the neutron and proton pairing transition temperatures 
T A for the different Sm isotopes. 

Fig. 8 The correlation between the transition temperature Ta and the ground state pairing 
gap A for different nuclei in the rare-earth region. 

Fig. 9 Temperature evolution of the hexadecapole moment (in units of barn 2 ) for the nuclei 
148 Sm and 150 Sm. 
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